On Positivity Preservation in Some Finite Element Methods for the Heat Equation
نویسنده
چکیده
We consider the initial boundary value problem for the homogeneous heat equation, with homogeneous Dirichlet boundary conditions. By the maximum principle the solution is nonnegative for positive time if the initial data are nonnegative. We study to what extent this property carries over to some piecewise linear finite element discretizations, namely the Standard Galerkin method, the Lumped Mass method, and the Finite Volume Element method. We address both spatially semidiscrete and fully discrete methods.
منابع مشابه
On Preservation of Positivity in Some Finite Element Methods for the Heat Equation
We consider the initial boundary value problem for the homogeneous heat equation, with homogeneous Dirichlet boundary conditions. By the maximum principle the solution is nonnegative for positive time if the initial data are nonnegative. We complement in a number of ways earlier studies of the possible extension of this fact to spatially semidiscrete and fully discrete piecewise linear finite e...
متن کاملNonstandard explicit third-order Runge-Kutta method with positivity property
When one solves differential equations, modeling physical phenomena, it is of great importance to take physical constraints into account. More precisely, numerical schemes have to be designed such that discrete solutions satisfy the same constraints as exact solutions. Based on general theory for positivity, with an explicit third-order Runge-Kutta method (we will refer to it as RK3 method) pos...
متن کاملMoving Mesh Non-standard Finite Difference Method for Non-linear Heat Transfer in a Thin Finite Rod
In this paper, a moving mesh technique and a non-standard finite difference method are combined, and a moving mesh non-standard finite difference (MMNSFD) method is developed to solve an initial boundary value problem involving a quartic nonlinearity that arises in heat transfer with thermal radiation. In this method, the moving spatial grid is obtained by a simple geometric adaptive algorithm ...
متن کاملFinite Element Methods for Convection Diffusion Equation
This paper deals with the finite element solution of the convection diffusion equation in one and two dimensions. Two main techniques are adopted and compared. The first one includes Petrov-Galerkin based on Lagrangian tensor product elements in conjunction with streamlined upwinding. The second approach represents Bubnov/Petrov-Galerkin schemes based on a new group of exponential elements. It ...
متن کاملPositivity-preserving nonstandard finite difference Schemes for simulation of advection-diffusion reaction equations
Systems in which reaction terms are coupled to diffusion and advection transports arise in a wide range of chemical engineering applications, physics, biology and environmental. In these cases, the components of the unknown can denote concentrations or population sizes which represent quantities and they need to remain positive. Classical finite difference schemes may produce numerical drawback...
متن کامل